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Abstract 



High precision cosmological observations in last decade suggest that about 70% of our 
universe's energy density is in so called "Dark Energy" (DE). Observations show that DE 
has negative effective pressure and therefore unhke conventional energy sources accelerates 
the cosmic expansion instead of decelerating it. DE is highly uniform and has become a 
dominant component only recently. 

The simplest candidate for DE is the time-independent cosmological constant A. Al- 
though successful in fitting available data, the cosmological constant model has a number 
of theoretical shortcomings and because of that alternative models of DE are considered. In 
one such scenario a cosmological scalar field that slowly rolls down its potential acts like a 
time-dependent cosmological constant. 

I have used different independent cosmological data sets to constrain the time dependence 
of DE's energy density in the framework of the slowly-rolling cosmological scalar field model. 
Present data favors a time-independent cosmological constant, but the time-dependent DE 
can not be ruled out at high confidence level. Ongoing and planned cosmological probes and 
surveys will provide more and better quality data over the next decade. When the new data 
sets are available we will be able to either detect the time dependence of DE or constrain it 
to a very small physically uninteresting value. 
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Chapter 1 
Introduction 

1.1 Short overview of relevant topics from theoretical 
cosmology 

Observations suggest that our Universe started evolving from an extremely hot and dense 
state 13 to 14 Gyr ago. Since light travels at finite speed, the distance to which we can 
see is fundamentally limited. The observable Universe today has a Hubble radius of about 
4000 Mpc, beyond which we can not see. Astronomical observations show that, on average, 
our Universe is homogeneous and isotropic on large scales. Observational evidence for 
isotropy includes the measurement of cosmic microwave background radiation anisotropies 
[1], shown in Fig. |1.1[ that is very isotropic. 

Large-scale structure observations confirm that the galaxy distribution becomes homo- 
geneous beyond about 100 Mpc, but even on the largest observed scales inhomogeneities in 
galaxy density are on the order of a few percent [2]. 

Homogeneous and isotropic four dimensional space-times are described by the metric 

ds^ = g^ydx^dx" = —dt^ + a{t)'jijdx^dx-' , (1-1) 

where t is cosmic time and jij is the metric of the three-space which can be fiat, spherical, 
or hyperbolic. a{t) is called a scale factor and describes how the spatial separation between 
two noninteracting test particles changes with time. In an expanding Universe the scale 
factor is an increasing function of cosmic time. 

1 



Figure 1.1: Map of cosmic microwave background anisotropics as measured by the Wilkin- 
son microwave anisotropy probe (WMAP). Red areas correspond to hotter spots on the sky, 
while blue areas correspond to colder than average spots. The amplitude of fluctuations 
relative to the background is less than 10~^. Source: http://map.gsfc.nasa.gov/. 



In general relativity (GR), the metric obeys Einstein's equations 



Rfiu + -^gtJivRx = -^T^y, (1.2) 



where G is Newton's gravitational constant, Rny is the Ricci tensor 



5rL d^\u 

P P'^ I pA pK pA pK (-1 q\ 



and the Christoffel symbols F^^^ are given by 



1 A« 99Kf^ , dg^^ dg^^ 



'^'^ 2^ I aa;'^ aa;/^ aa;'^ ' ■ ^ ' ^ 



T^i, in Eq. (1.2) is the energy-momentum tensor of the matter components that fill space. 
The Ricci tensor satisfies the Bianchi identity -R(^.^ — 1/2R^.^ = 0, where the semicolon 
denotes a covariant derivative, and consequently the energy-momentum tensor satisfies the 
energy conservation equation 

= 0. (1.5) 

Three-dimensional homogeneous and isotropic space can be fiat, positively curved (sur- 
face of a 3-dimensional "sphere"), or negatively curved (surface of a 3-dimensional "saddle"); 
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two dimensional analogs are shown in Fig. 



1.2 
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Figure 1.2: Examples of 2 dimensional flat, closed, and open spaces. Source: 
http://www. Icsd. gov. hk/. 



The metric jij of homogeneous and isotropic space can be written in Cartesian coordi- 
nates as 



'jijdx'^dx^ 



^ 1 J doc doc"^ 
(1 + iKp^) 



2 ' 



1.6) 



and in spherical coordinates as 



7.^. = dr^ + x\r){d9^ + sm\9)df 



;i-7) 



Here 6ij is the Kroneker delta function and = ^jjX^X"' . K describes the geometry 
of three-space. For K < the space is negatively curved (open), for > the space is 
positively curved (closed), and for if = the space is flat. The variable x is given by 



sin \/Kr 



/K 



X[r) 



if is: > 0, 

if K = 0, 
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^ sinh ^/^r if K < 0. 



The numerical value of K depends on the normalization of the scale factor. In open and 
closed Universes a is usually normalized in such a way that K = ±1. In the spatially-flat 
Universe the scale factor is often normalized so that it's current value is equal to 1. 



If the metric is given by Eqs. (1.1) and (1.7) than Eq. (1.2) forces the energy-momentum 
tensor to have the form 



( p{t)goo \ 






;i.9) 



where p is the mass energy density of matter and P is pressure. 

For barotropic fluids pressure depends only on the local value of energy density, P = 
P{p). Very often the equation of state P = up is used, where u is called the equation of 
state parameter. For an ultrarelativistic gas (radiation) uj^ = 1/3, for nonrelativistic matter 
cUm = 0, and for a cosmological constant u/^ = —1. 



Energy conservation, Eq. (1.5), in the expanding Universe gives 



p = -3(1 + u)p- 



;i.io) 



For a general barotropic fluid Eq. (1.10) can be solved and yields 

'ao\ 3(1+'-') 



P = Po 



a / 



;i.ii) 



;i.i2) 



which for uj = —1 becomes 

p = po {uj = -1), 

where po is the current energy density of the fluid. 

The combination H = a/a is known as the Hubble parameter, where an overdot denotes 
a derivative with respect to cosmic time t. It shows how fast the Universe is expanding. Its 
current value Hq is called the Hubble constant and is often quoted in terms of a dimensionless 
parameter h = H/{100 km s~^ Mpc"^). 



In the Universe described by Eqs. (1.1) and (1.7) and dominated by a single barotropic 



fluid, Eq. (1.2) reduce to two independent equations for the scale factor 
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Equation (1.10) follows from Eqs. (1.13) and (1.14). In the Universe dominated by a single 



barotropic fluid, Eqs. (1.13) and (1.14) give 

/ ^ \ 2/3(1+0.) 



(c^-l). 



;i.l5) 



Equation (1.13) shows that if p = pcr = 3H^/8ttG the curvature is zero. This time-dependent 
energy density that makes the Universe spatially-flat is called the critical density. Current 
energy densities of matter components are often given in terms of dimensionless density 
parameter, Qi = Pi/pcr- 

In a curved space with radiation, nonrelativistic matter, and cosmological constant 



Eq. (1.13) can be rewritten as 



;i.i6) 



where Q^a, ^r, ^a, and are respectively the density parameters of nonrelativistic matter, 
radiation, cosmological constant, and curvature, and the scale factor is normalized so that 
ao = 1. 

1.2 Distance and time measurements in cosmology 

Distances can be measured directly, using the parallax method, only to astronomical objects 
that are close to us (less than 100 pc). For objects that are farther away indirect methods 
of distance determination must be used. 

When a distant object emits a photon with a wavelength Agm at time tem, since light 
travels at finite speed, we will detect the photon at a later time tnow As the Universe 
expands, the photon will "redshift" and we will detect it at a longer wavelength 



A - A 



1.17) 



where aem and ao are respectively scale factors at the moment of emission and detection 
(now) of the photon. Redshift of an astronomical object is defined as 

ao 



^obj 



- 1 



C^obj 
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1.18) 



where Oobj is the value of the scale factor at the moment when the light that reached us 
from that body was emitted. 

If the cosmological expansion history a{t) is known, redshift can be uniquely related to 
the distance to the object and the time of light's emission. Distances and times in cosmology 
are often quoted in terms of redshifts. Two most frequently used distance measures in 
cosmology are the luminosity distance and the angular distance. 

1.2.1 Measuring time intervals 

Light travels on null geodesies cis^ = 0. When we detect a photon that was emitted at 
redshift z it has traveled for time AT = D/c, where D is the distance and c is speed of 
light. Let's set up the coordinate system so that the photon travels in the radial direction. 



Then from Eqs. (1.1) and (1.7) the total elapsed time can be expressed as 



rtiz) j-a{z) 

AT = dt= ax{a)da. (1.19) 

Jt{0) Jao 



Using Eq. (1.16) and the identities a = ao/(l + z) and da = —dzao/{l + z) , we can rewrite 



Eq. (1.19) as 



AT = t{z) - to = r (1-20) 



In a given cosmological model Eq. (1.20) uniquely relates redshift to time. 



1.2.2 Luminosity distance 

Let's say we have a source of known luminosity (power) L at a redshift z. The source emits 
energy Ldt in a time interval dt and this energy redshifts hj a/a^ = 1/(1 + z) before it 
reaches us. The energy is distributed over a spherical surface with radius aox(^o — t), where 
x{to — t) is the distance that a photon travels before detection. The flux that we measure 
in unit time interval dto = {ao/a)dt is given by 



47ra^x'(to-t)' 



;i.21) 



Luminosity distance is defined as 



AtiF 



a[z) 



(1.22) 



By substituting to ~ from Eq. (1.20), the luminosity distance can be rewritten as 

dz 



di 



Ho Jo ^/Vt^{l + zf + a(l + zY + ^K + ^k{1 + 



(1.23) 



Luminosities of some astronomical objects can be determined by indirect measurements. 
Supernovae Type la (SNela), for example, have very small dispersion in total luminosity. 
Gamma-ray bursts (GRBs) could be well calibrated in the future with their luminosity 
determined to high accuracy. We can measure the flux of these and other "standardizable 



candles" at different redshifts to determine cosmological parameters from Eq. (1.23) 



1.2.3 Angular diameter distance 

Let's say an astronomical body at a distance d from us has a physical size D, then the 
observed angle subtended by that body is 9 = D/d. The angular diameter distance is 
defined as 



d. 



D 
J' 



(1.24) 



If the object is at redshift z it will span the angle 6 = a{z)x{z). Using arguments of 



Section 1.2.2 we can express the angular diameter distance through cosmological parameters 
as 



d 



1 



A 



-X 



dz 



1.25) 



1 + z \aoHoJo + z)^ + ^^r(l + z)^ + ^^a + ^^k(1 + z)^ 

Similarly to standard candles, "standard rulers", such as FRII radio galaxies, have small 



dispersion in physical size. Equation (1.25) then can be used to determine cosmological 



parameters from the redshifts and apparent angular sizes of those standard rulers. 
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1.3 Short overview of relevant topics from observa- 
tional cosmology 



In last decade the quality and quantity of astronomical observations relevant to cosmology 
have greatly improved. SNela can map the expansion of the Universe up to the redshift 



of 1.7 (see Fig. 1.3). While the power spectra of CMB anisotropics (see Fig. 1.4) and the 



galaxy distribution (see Fig. 1.5) have been measured to a high precision. 




Figure 1.3: Hubble diagram of SNela detected by ESA-WFI. Black dots represent SNela 
with measured absolute magnitude at different redshifts. Source: http://www.sci.esa.int/. 



The combination of SNela absolute magnitude versus redshift data [3], E] , measurements 
of CMB anisotropy [5], cluster gas mass fraction versus redshift data |6j, measurements of 
the baryon acoustic peak [71 [8], and other large-scale structure measurements [9] can be 
used to determine cosmological parameters with high precision. 

Currently available cosmological data strongly indicate that: 

• The Universe is highly isotropic and homogeneous. 

• The Universe is expanding and distance between gravitationally unbound objects are 
increasing. 
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Figure 1.4: CMB temperature anisotropy power spectrum. Black dots with error bars denote 
different measurements, red line shows predictions of the best-fit theoretical model. Source: 
http: //cosmos, astro. uiuc. edu/. 



• The cosmic expansion started accelerating at a redshift near z = 0.5, if the dark energy 
is a cosmological constant. 

• The spatial curvature of the Universe is very close to zero, if the dark energy is 
cosmological constant. 

The accelerated expansion of the Universe can be understood if we assume that around 
75% of the current cosmological energy budget is made of dark energy (DE). Parameterized 
as a fluid, DE has a negative equation of state parameter < ^1/3 and is spatially 
uniform. 

About the other matter components we know that: 

• The Universe is filled with CMB photons. The CMB has an accurately Plankian 
spectrum with a temperature of T = 2.73 K. Radiation today makes a negligible 
contribution to the total energy density. 

• Baryons make up less than 5% of critical density. About 75% of baryonic matter is 
Hydrogen, 25% Helium. Contribution of heavier elements to the total baryonic energy 
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Figure 1.5: Power spectrum of giant red galaxies as measured by Sloan Digital Sky Survey. 
Black dots with error bars denote measurements, different colored lines show theoretical 
predictions. The baryon acoustic peak is zoomed in on in the upper right corner insert. 
Source: http://www.sdss3.org/. 

density is small. The entropy of the Universe is high, around 10^ photons per baryon. 

• Nonrelativistic, cold dark matter (CDM) contributes around 25% of the critical den- 
sity. CDM is not baryonic and does not significantly emit, absorb, or scatter light. 

A joint analysis of current CMB anisotropy measurements, SNela data, and BAO peak 
measurements results in the following la ranges for key cosmological parameters [TO] : 

• h = 0.71 ±0.013. 

• f^b = 0.0462 ±0.0015. 

• ^DM = 0.233 ±0.013. 

• = 0.721 ±0.15. 

• Age of the Universe = 13.73 ± 0.12 Gyr. 

This chapter provides only a concise overview of topics relevant to the rest of this thesis. 
Detailed, in depth treatments of theoretical and observational cosmology can be found in 
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various textbooks [TH4T5] . 



Chapter 2 
Dark Energy 



2.1 Observational evidence for dark energy 

The first direct evidence of accelerated expansion was provided by SNela absolute magnitude 
versus redsliift data [151 HT] . Accelerated expansion is very difficult to explain in GR with 



only conventional matter. Equation (1.14) implies that the expansion can be accelerating 
only if the Universe is dominated by a component with negative pressure P < — p/3. The 
initial SNela measurements were followed by ground-based and space surveys, such as the 
Hubble Space Telescope (HST) survey [IH], the supernova legacy survey (SNLS) [IH], and 
the ESSENCE survey [20]. Current SNela data extends to z = 1.7 and provides more than 
5(T evidence for accelerated expansion. 

Positions and amplitudes of acoustic peaks in the CMB anisotropy show that the Uni- 
verse is very close to spatially flat if DE does not evolve in time. When CMB anisotropy 
results are combined with an independent measurement of the Hubble constant they in- 
dicate that nonrelativistic matter makes up only about a quarter of the critical density, 
while the energy density of radiation is negligible. This implies that about three quarters of 
our Universe's energy budget is made of something different than ordinary or dark matter 
and radiation. This extra component must be uniform on cosmological scales and became 
dominant recently, in order not to adversely effect large-scale structure formation and big 
bang nucleosynthesis. 
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The presence of DE also effects CMB anisotropies directly through the integrated Sachs- 
Wolf (ISW) effect. In a Universe dominated by DE gravitational wells decay on large dis- 
tance scales. This affects the large-angle CMB spectrum, since photons that travel through 
decaying gravitational potentials gain energy. If the Universe is dominated by DE there 
should be cross-correlation between CMB anisotropies and low-redshift matter distribution. 
Several groups reported detection of the ISW effect in cross-correlations of the CMB with 
galaxy and radio-source catalogs [2TH23] . 

These results, combined with the BAO peak measurements, LSS tests, and galaxy cluster 
data, provide fairly overwhelming evidence that about 70% of the current Universe's budget 
is in spatially uniform DE with negative pressure. 

2.2 ACDM model 

All currently available cosmological data is quite consistent with a Universe that is spatially- 
flat, is dominated by a cosmological constant A with Q\ ~ 0.7, with the rest of the energy 
density being in nonrelativistic cold dark matter with ~ 0.25 and nonrelativistic baryonic 



matter with fib ~ 0.05. A cosmological constant has uja = —1 and according to Eq. (1.10) 
it's energy density does not change in time. This model is referred to as the "standard 
model" or the spatially-fiat ACDM model. 

In the spatially-fiat ACDM model the background expansion of the Universe at late 
times is described by 

H = Ho^/n^ii + z)3 + nA. (2.1) 

Although the ACDM model is a good fit, the data is not yet tightly constraining and does 
not yet require a constant A. To account for the possibility of a time-varying cosmological 
constant the equation of state of DE is sometimes modeled as P = uj^p, where is a negative 
parameter. This, so called, XCDM parametrization lacks physical motivation. XCDM is 
usually used as an ansatz in data analysis for the purposes of quantifying time-dependence 
of DE. 
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2.3 Problems of ACDM 



Although successful in fitting most data, the ACDM model has a number of theoretical 
shortcomings. A cosmological constant is difficult to motivate from fundamental physics. 
The most plausible candidate for a cosmological constant is vacuum energy. 

Since all Lorentz observers should see the same vacuum, it's energy-momentum tensor 
must be proportional to the metric tensor TJ^^ = AgJ^^. The vacuum therefore has the 
required equation of state P™'^ = p™'^. It's very difficult however, to explain the small value 
of DE energy density based on the vacuum energy model. 

In the quantum field theory the zero point energy of the vacuum is divergent 



f oo 

p^'^oc / k^dk. (2.2) 
Jo 

Quantum field theory is expected to break down at the Planck scale of around 10^^ Gev. 



If we use this as upper integration limit in Eq. (2.2) (instead of oo) we get a huge number 
that exceeds the observed dark energy density by 120 orders of magnitude. 

In a supersymmetric model every boson has a fermion of equal mass as a supersymmet- 
ric partner and the vacuum energies of these partners cancel. Supersymmetry (SUSY), if 
existent, is believed to be broken at an energy of roughly 1 TeV or so. If we cut off the upper 



integration limit in Eq. (2.2) at the energy of SUSY breaking we will still get a difference 
of around 60 orders of magnitude. This discrepancy between the small measured value of 
cosmological constant and the much larger theoretically "expected" values of vacuum energy 
is known as the "smallness" problem |24j . 

One possible explanation of the "smallness" problem is based on anthropic arguments. In 
string theory, multiple vacuum states with all possible values of vacuum energy are possible. 
Different causally disconnected patches of the Universe spontaneously choose vacuum states 
that are independent of each other. If the Universe is infinite there will always be parts of it 
that have a given value, no matter how unlikely, of the vacuum energy and we just happen 
to live in one of those regions with a very small value of vacuum energy density [251 126] . 
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Another interesting fact that's difficult to explain in the ACDM model is that today both 
nonrelativistic matter and DE have comparable energy densities. This is surprising since the 
matter and DE components scale with redshift differently. For radiation it is ~ (1 + z)^, 
for CDM and baryons ~ (1 + z)^, and for a cosmological constant Q/<^ ~ const. 

At the beginning of cosmic evolution the Universe was radiation dominated, today ra- 
diation contributes less than 1% of the total energy density. The contribution of DE was 
negligible in the past, it has become a dominant component only recently, and in the future 



will be the only component driving cosmic expansion as shown in Fig. 2.1 
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Figure 2.1: The time evolution of energy densities of cold dark matter (green line), radia- 
tion (red line), and cosmological constant (black line). 



There is only a short period of time when the energy densities of matter and cosmological 
constant are comparable. It is unclear why we happen to live in this narrow window of time. 
This is called the "coincidence" problem and like the "smallness" problem can be addressed 
by the anthropic arguments. 

Besides the two problems mentioned above there are other observational facts that ap- 
pear to conflict with the predictions of the ACDM model, at possibly more then 2a confi- 
dence level. These are 

• High redshift SNela data are consistent with spatially-fiat ACDM. It favors however 
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models with ujbe < — 1- Initially, this was thought to be a statistical fluke that would 
go away as more data accumulated. This discrepancy however still persists even as 
larger data sets become available. The discrepancy is caused by high {z > 1) redshift 
SNela which are systematically brighter than what we would expect in the ACDM 
model [3]. This could be due to an unknown systematic effects, possibly associated 
with high-redshift SNela evolution, or a statistical effect that will go away with more 
data. If the discrepancy persists, it would mean that the Universe in the past was 
decelerating faster than the ACDM model predicts. 

• Large-scale velocity flows have amplitude of 400 km/s, larger than what is expected 
in a ACDM model [27H3T] . Velocity flows extend to z = 0.2 and could be as large 
as 1600 km/s. In ACDM the probability of having velocity flows with such a large 
amplitude is less then 1%. The explanation for this could be that it is just a big 
statistical fluctuation or there could be some physical reason, such as time-dependent 
Newton's constant, presence of noninflationary perturbations, or a giant void at a 
distance of few Gpc. 

• Cosmological simulations based on the ACDM model predict that large voids should be 
filled with many dwarf dark matter halos. This turns out to be true for very large voids 
(larger than 10 Mpc). Smaller voids however are observed to be surprisingly empty 
of dark matter halos [32]. For example, based on ACDM we would expect to observe 
on average 10 dwarf galaxies in our local void, but there are none. Possible resolution 
of this problem could be related to the incompleteness of observational sample, or 
incorrect bias model that fails to account for specific environmental properties. 

At the moment these discrepancies between observations and theoretical predictions do 
not have big weight in the overall picture. The ACDM model is in general a good statistical 
fit to the combined data. These inconsistencies might go away when new high-quality data 
become available. 
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2.4 Alternative dark energy models 

At the moment, ACDM is a good fit to available cosmological data, however, because of 



the observational issues mentioned in Sec. 2.3 and its inability to naturally account for the 



small observed energy density of DE, a number of alternative models have been proposed 
over the years. 

One approach is to not introduce a new exotic DE component, but to modify the equa- 



tions of GR and consequently Eq. (1.2), [33H37] . In this set of models gravity is usually 
weaker on larger scales and the accelerated expansion is just a gravitational effect. At very 
early times modified gravity should act like GR predicts, in order not to conflict with pre- 
dictions of CMB anisotropy and BBN. In f{R) and scalar-tensor theories, the equivalents 



of Eq. (1.2) include up to forth order derivatives of the metric tensor. Gravity is not a 
spin-2 transverse and traceless field anymore. It is also carried by spin-0 (scalar) and spin-1 
(vector) particles. 

Another possibility is that the Universe is not as homogeneous as we thought and the 
accelerated expansion is an effect of averaging nonlinear spatial inhomogeneities over cos- 
mological scales [381 - HTj . It is not yet clear if the backreaction of structure formation on 
the background evolution can be strong enough to mimic accelerated expansion and at the 
same time be consistent with all other available data. 

A third option is to introduce a DE fluid that is close to spatially uniform but unlike a 
cosmological constant slowly varies in time. This is usually achieved by introducing a scalar 
field that is slowly rolling down it's almost flat potential [l2l |33]. In the minimal model, 
the scalar field is coupled to the rest of the matter only gravitationally, but other couplings 
have also been considered jUl - liT] . 

In Chapter |3] I will elaborate more on the scalar field dark energy models. More detailed 
discussion of DE theory and the observational situation can be found in recent reviews 
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Chapter 3 

Cosmological Scalar Field 



3.1 Scalar field in the expanding universe 

Let's assume that a self-interacting scalar field (j) minimally coupled to gravity is present on 
cosmological scales. This scalar field is described by the lagrangian density 

c = ]^d,ct>d^<p-v{<p), (3.1) 

where V^(0) is the potential energy density. The energy momentum tensor of this scalar 
field is given by 

d(t> d(f) „ fl d(f) (90 ,A , , 

If the scalar field is close to spatially uniform on cosmological scales we can neglect its 
spatial derivatives (p^i compared to time derivatives and the individual components of the 
homogeneous part of the energy momentum tensor can be written as 

(3.3) 
(3.4) 
(3.5) 
(3.6) 



rjiO 

i 


= 0, 




rpi 
j 


= 




rpO 

-'o 


= p -- 




rpi 

i 


= p 





By looking at Eqs. (3.5) and (3.6) we can see that the equation of state parameter 



.^€^^ (3.7) 
02 + 2V{ct)) 
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can become negative if the scalar field evolves slowly in time. The equation of state pa- 
rameter is always greater than —1 and in the limit of the scalar field "at rest" (0 = 0) we 
effectively get a cosmological constant with u = —1. 

The time evolution of the scalar field (and consequently the time evolution of the equation 
of state parameter) is determined by the Klein-Gordon equation. For the metric given by 



Eq. (1.1) the Klein-Gordon equation is 

^ + 3^^+?^=0. (3.8) 

a 0(j) 

If the scalar field energy density makes nonnegligible contribution to the total energy budget 



of the Universe it will affect the evolution of the scale factor and Eq. (1.16) should be 
rewritten as 

H = HJ^ + ^ + Q,{a) + ^, (3.9) 

where Qfp is now the time- dependent energy density parameter of the scalar field. 

To make specific predictions about how the scalar field affects the Universe's evolution 
we have to specify the form of the potential V{(j)). Since the underlying physics is not 
known at the moment, it's impossible to pick a specific functional form for 1^(0) based on 
the considerations of fundamental physics only. 

The only fundamental scalar field in the standard model of particle physics is the Higgs 
field. Rest mass of the Higgs boson is larger than about 100 GeV and at low energies it 
decays into W bosons, H W~^W~ , therefore Higgs bosons can not be present in large 
quantities in the late universe. In the theories that go beyond the standard model, such as 
string theory, minimal supersymmetric theory, etc., fundamental scalar fields arise naturally 
and we have a large number of candidates. The particle physics at very high energies is not 
well constrained experimentally yet and different theoretical scenarios are possible; because 
of this uncertainty its impossible to motivate one choice of potential based on fundamental 
considerations and a number of different models have been proposed over the years. An 



incomplete list of possible scalar field potentials is shown in Table 3.1 
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Potential V{(f)) 


Reference 


0"°, a > 


US] 


exp(— A0) 


m 




[52] 


exp(-A0)/0", a > 




exp(A/0) — 1 


m 


(cosh(A0) - > 




((0-Ci)^ + C2)exp(-/c0) 


[56] 



Table 3.1: List of proposed self-interaction potentials for the cosmological scalar field. 



In this work I will concentrate on the inverse-power law potential first proposed in 
Ref. ^42j. In this model the potential is assumed to be inversely proportional to a power of 
the scalar field, V oc 1/0", where a is a positive constant. 

3.2 Inverse-power law potential 

If we take scalar fields self-interaction potential to be inversely proportional to a power of 



V^(0) = ^0-", (3.10) 



where k and a are nonnegative parameters, Eqs. (3.9) and (3.8) can be rewritten as 



SirG 



a J 3 



(Pm + P</,), (3.11) 



■■ a ■ K a , ^ 

where we have set spatial curvature to zero and neglected radiation, and the energy density 
of the scalar field is given by 

P*^^^^+^). (3.13) 

Parameter a describes the steepness of the scalar field potential. Larger values of a cor- 
respond to faster evolution of the scalar field and vice versa, a = corresponds to the 
time- independent cosmological constant. Parameter k sets the mass scale M of the scalar 
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particle, 



/ K \ -("+4) 

M-{^) . (3.14) 



In the distant past, when the energy density of the scalar field is much less than the 



energy density of other components <^ Pcdm , Eq- (3.12) accepts a solution 

(j) oc a3(i+-cd„.)/("+2)^ (3^15) 

and the equation of state parameter of the scalar field is 

= — (3.16) 

For positive a, redshifts less rapidly than the dominant component of energy and at some 
point eventually starts to dominate. The redshift when cosmic acceleration starts is given 
by 



As we approach zs-, the fractional energy density of the scalar field increases and the ap- 



z^=^ -1. (3.17) 



proximation p^ <^ Pcdm breaks down. Evolution of the scalar field does not follow Eq. (3.15) 



anymore and to solve for (j)[a) now the set of coupled differential equations (3.11) and (3.12) 
must be solved numerically. 



A nice property of Eq. (3.15) is that it is an attractor. A wide range of initial conditions 



on (j) result in solutions that converge to the solution of Eq. (3.15) at some point. Figure 3.1 
shows the evolution of the scalar field for different initial conditions. Although initial con- 
ditions on p^ span some 70 orders of magnitude, the late-time evolution of the scalar field 
is similar in all intermediate cases. 

In the scalar field scenario the "smallness" and "coincidence" problems mentioned in 
Sec. [2] are partially solved. The scalar field can start from a very high energy state, roll 
down the inverse power law potential, and naturally lead to a that evolves to the very 



small measured value. Figure 3.1 shows that the scalar field "tracks" the evolution of 



the dominant component once it joins the solution of Eq. (3.15). The energy densities 
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Figure 3.1: Evolution of energy densities of nonrelativistic matter (green line), radiation 
(red line), and scalar field with different initial energies (black lines). Scalar fields that start 
off from different energies have very similar late-time evolution. 

of nonrelativistic matter and scalar field are comparable for a longer period of time, and 
current time becomes less special. 

The scalar field in this model does not cluster on scales much smaller than the Hubble 
radius 1/Hq as required by CMB and large-scale structure data. 



In our analysis we will assume that the only two constituents affecting the late-time evolu- 
tion of the Universe are nonrelativistic cold dark matter and baryons with energy density 
parameter fi^ and a scalar field with inverse power law self-interaction potential V cx 
We will also assume that the Universe is spatially-fiat which implies + = 1- 

In this model the background expansion is fully described by two parameters, f2m and 
a. The evolution of the scale factor and scalar field are governed by 



3.3 (/)CDM model of dark energy 




(3.18) 



+ 3-0 + 



a 



dV{4>) 
d(j) 



0. 



(3.19) 
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For a = the (pCDM model reduces to the spatially-flat ACDM with Qa = Q,^. For 
nonzero a predictions of 0CDM and ACDM differ, for example the scalar field model predicts 
lower values of luminosity and angular diameter distances for the same value of Hubble 
parameter and nonrelativistic matter density. An example of a C++ function that solves 
Eqs. (3.18) and (3.19) is presented in Appendix [A} 

Different cosmological tests have been used in the past to constrain the parameters of the 
0CDM model, including counts of galaxies [57], angular sizes of radio sources and quazars 
[SHI EH], galaxy cluster gas mass fraction data [HD], strong gravitational lensing pT] and 
SNela absolute magnitude versus redshift data [62146^ . In next chapters I will discuss the 
constraints on 0CDM model from some current data sets. 

For an in depth review of scalar field dark energy see, e.g., Refs. [121 113|, [651 166] . 
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Chapter 4 

Data Analysis Techniques 



We are interested in observable quantities X°^^ measured at redshifts Zi (or in redshift bins 
of width Azi). X°^^ could be, e.g., luminosity distance, angular diameter distance, Hubble 
parameter, or quantities that implicitly depend on them. First I will assume that evolution 



values of parameters a and ^l^. For given values of model parameters we can compute the 
theoretical expectations X^^ for the observed quantities. We will find the best-fit parameters 
a* and Q^, for which the theoretical predictions X^^ are the closest match to the results 
of observations X°^^ by some measure. We will also construct the confidence level intervals 
that are likely to cover the true values of parameters with a specified probability. 

We will compare the predictions of 0CDM model to the ones of ACDM and see which 
ones fit the observations best. In the end we will see if data favors one model over the other, 
and quantify the degree of discrepancy between the two models. 



4.1 likelihood functions 

Let's say we have N measurements X°^^ at redshifts Zi with measurement errors di, and the 
0CDM model predicts Xl^{a,flrn)- The function. 



of the Universe is described by the 0CDM model, as presented in Chapter |3} with unknown 




(4.1) 
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quantifies the discrepancy between theoretical predictions and observations. Small values of 
indicate a good fit. The parameters a* and that minimize called the best-fit 

parameters. The likelihood function, 

C{a, nj = exp {-x^{a, QJ/2) , (4.2) 

has a local maximum Cml at {a*, ^^). If the measurements are independent and Gaussian 
distributed with mean X"^^ and variance (Ji, then the best-fit values of parameters are 
unbiased estimators of their true values. Values of parameters that result in high values of 
the likelihood function are more likely to be the true parameters. 

1, 2, and 3a confidence level contours are defined as the pairs of parameters {a, flra) for 
which — 2 ln£(«, f2m) is less by 2.30, 6.18, and 11.83 respectively compared to Cml- If the 
likelihood function had an exact two dimensional Gaussian profile, 

jC oc e-(<^i(°-°*)'+'^2(n^-nm)2)/(2a2)^ ^^ g-j 

Ncr contours would correspond to the range of parameters {a* ± A'"cr, f)^ ± Na). In that 
case the 3a contour would enclose 99.73% of integrated hkehhood and the parameter values 
outside that range would be 0.27% or less hkely to be the true parameters. This is not true 
for an arbitrary £, however, when it is a sharply peaked function, parameter values outside 
3a are usually less than 1% likely. 

4.2 Nuisance parameters 

Very often X^^{zi) will depend on parameters other than a and fl^. These "nuisance" 
parameters u can be determined from the data but are of no interest to us. Examples of 
nuisance parameters in dark energy parameter estimation are the Hubble constant Hq and 
baryonic energy density parameter f2b- In this case the likelihood function C{a,u!m,u) will 
also depend on i^*. If we have some idea from other observations about what the values of 
u are, we can summarize this information in the form of probability distribution function 
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(p.d.f.) P{i^). P{i^) is normalized to one and is usually taken to be a Gaussian peaked at 
the best-fit value of u* determined from other independent data sets. If we do not want 
this implicit use of information from other data sets, we can use the conservative flat prior 
P{h') = l/(i^2 — J^i) over the range (i/i,i/2). We then integrate the likelihood function, 

£'(«, n^) = J C{a, n^n, u)P{u)du. (4.4) 

and determine the best fit values and confidence level contours from C as in Sec. 14. 11 

Two most frequently used priors on the Hubble constant are the measurement of Hub- 
ble Space Telescope (HST) team |67] Hq = (72 ± 8) km s^^ Mpc~^ and the result from the 
Wilkinson Microwave Anisotropy Mission (WMAP) 5-year data [HH] -f^o = (73±3) km s^'^ Mpc~^. 
The WMAP measurement has less dispersion but is model dependent, unlike the HST mea- 
surement which measured Hq directly based on the recessional speeds of distant objects. 
The most conservative option is to integrate over Hq with a flat, noninformative prior. The 
baryonic matter energy density is also well constrained by a number of independent tests. 
Big bang nucleosynthesis gives the value i7b^^ = 0.0205 ± 0.018 [69j and the WMAP data 
results in n^h^ = 0.0223 ± 0.0008. 

4.3 Constraints on individual cosmological parameters 

One-dimensional confidence level contours and best-fit values can be computed for individual 
cosmological parameters a and Q^- We take the two-dimensional likelihood function from 



Eq. (4.4) and integrate it with respect to the other parameter with a flat prior. 

C^ = J C'ia,n^)dn^, (4.5) 

C,n = j C'{a,Q^)da. (4.6) 

(4.7) 

For each parameter we flnd the best-flt value that maximizes the corresponding one- 
dimensional likelihood function. We deflne 1, 2, and 3a confldence level intervals {ai,a2) 
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and (r2mi,f^m2) so that 

s = j Cada / / Cada, (4.8) 

J ai / J all a 

s= / Cn^^dn^/ / Cn^^n^, (4.9) 

J flml I J all 

and s equals to 68.27, 95.45, and 99.73 respectively. We will choose lower and upper limits 
on parameter values so that the likelihood is higher everywhere inside the interval than 
outside. 

Parameter a should always be positive. Since the best-fit values are usually close to the 



a = line, ol\ in Eq. (4.8) will often be equal to zero. 
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Chapter 5 



Constraints on ^CDM from Hubble 
Parameter Measurements 



This chapter is based on Refs. [701 - I72] . 

5.1 Hubble parameter measurements 

The Hubble parameter in spatially fiat 0CDM is given by 



H = H^^n^ + %{a) (5.1) 

and is a function of redshift z and cosmological parameters fim and a. In practice, it can be 
measured indirectly at different redshifts and thus used to constrain cosmological parameters 
of interest. 

First of all let's note that since a ~ (1 + z)^^, the Hubble parameter can be rewritten as 

^=- = -T^?- (5-2) 



Equation 5^ shows that the Hubble parameter can be determined from the measurement 
of dz/dt. dz/dt itself can be determined from a survey that measures redshifts of galaxies 
with high precision. The basic idea is that if we have two galaxies very close in redshift 
space and we can determine their ages with good accuracy, then dz/dt can be approximated 
by reasonably well with Az/AT, where Az is the difference in their redshift and AT is the 
difference in their age. The most difficult part is to accurately determine ages of galaxies. 
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In passively evolving galaxies the star formation rate is low and their ages can be deter- 
mined with better accuracy. These galaxies are dominated by the light from old red main 
sequence stars. The evolution of those stars is a well understood process and metallicity 
can be simulated on computers. The spectra of galaxies depends mainly on metallicity and 
age. A computed spectrum can be compared to the observed one and the galaxy age can be 
determined. Stars in galaxies, of course, do not have the same metallicities and this spread 
in metallicity introduces a statistical error of about 0.1 Gyr. For a detailed description of 
this and other relevant issues and methods see Refs. 



5.2 Constraints from Hubble parameter data 

Reference [71] used Gemini Deep Deep Survey [75] and archival data [75H5U] to estimate ages 
of 32 passively evolving galaxies. From these age measurements they were able to compute 
the differential age dz/dt at 9 different redshifts in the redshift range 0.09 < z < 1.75. The 
age measurements are shown in Appendix [Bj while the Hubble parameter estimates from 
Ref. [^ are shown in Appendix [C] These H{z) measurements are independent and the 
errors are Gaussian distributed. The function is 

X =2^ -2 (5.3) 

i=i « 



where i^th is the theoretical value of Hubble parameter computed from Eq. (5.1), and 
Hohs{zi) and ai are the observed values and one sigma uncertainties. In our analysis we 
assume 12% uncertainties on age measurements in Appendix [B] [8T] . 

The Hubble constant (the Hubble parameter now, at zero redshift) is a nuisance param- 
eter and can be integrated over to get a two dimensional likelihood function of cosmological 
parameters 



CiSlr^.a) = J exp{-x\n^,a,Ho)/2)P{Ho)dHo (5.4) 

where P{Ho) is a prior probability distribution function for Hq. We integrate over Hq with 
a Gaussian WMAP prior with Hq = (73 ± 3) km s^^ Mpc'^ 
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The 1, 2, and 3a contours on 0CDM are shown in Fig. 5A_ Data constrain the non- 
relativistic matter density parameter to be less than 0.4 at about 3a confidence level. The 
likelihood in the direction of a is almost fiat. The one dimensional likelihood functions on 



individual dark energy parameters are shown in Figure 5.2 a is poorly constrained while the 



nonrelativistic matter density parameter is constrained to be in the range 0.09 < Q^i < 0.26 
at la. 




0.45 



Figure 5.1: 1, 2, and 3a confidence level contours on (pCDM model parameters from Hubble 
parameter versus redshift data. 
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Figure 5.2: One dimensional likelihood functions of individual cosmological parameters 
from Hubble parameter versus redshift data. The maximum likelihood values are normalized 
to one. 

5.3 Lookback time 

The data in Appendix [B] can be used to directly constrain dark energy parameters using the 
lookback time versus redshift relation [ISIIS^HS^ . Lookback time at redshift z is defined as 
the difference between the current age of the Universe and it's age at redshift z. In 0CDM 
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the lookback time is given by 

t,^(z) = ^ r , =■ (5.5) 

For an object at redshift z it's age is the difference between the Universe's age at z and the 
one when the object was formed at redshift Zf. 

1 /"^f dz 
-no (1 + z)A/nm(l + + ^^0(2) 



From Eqs. (5.5) and (5.6) the observed lookback time is defined as 

C(z)=C-WW-^, (5.7) 

where t^^^ is observed total age of the Universe, t{z) is the objects age at redshift z, and 
r = Iq^^ — tib^f is a "delay factor" that encodes our ignorance of Z{. The "delay factor" 
T is different for each object in the sample. From CMB experiments the total age of the 
Universe is estimated to be tg^*^ = (13.7 ± 0.2) Gyr [68j. 
We construct the function 

X [a,\lyn, no,T) = 2 ^ 2 ' y^-^) 

i=l ^to 

and compute the likelihood function 

C{a, fim, t) oc exp(-x^(a, ^m, r)/2). (5.9) 
We account for our ignorance of r by marginalizing over it by integrating the likelihood 



function in Eq. (5.9) analytically over r with a fiat prior. 
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C\a,n^,Ho) = I Cdr = erfc (^^) exp {^^^^ ] ■ (5-10) 
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where 



i=l 



B = f;Mf!lz<!M, (5.12) 



i=l 
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C = X;4 (5.13) 



= (5.14) 

and erfc(x) is the complementary error function. Finally we integrate over the Hubble 
constant with a Gaussian WMAP prior with Ho = (73 ± 3) km s'^ Mpc^^ to get the 
constraints on (pCDM model parameters a and Q^. 

5.4 Constraints from lookback time versus redshift data 



The constraints on 0CDM from the lookback time method are shown in Fig. 5.3 

Lookback time versus redshift data constrains the nonrelativistic matter density to be 
less than fim = 0.35 at about 3a confidence level. The a parameter is not well constrained 
and values as large as a = 10 are allowed at la confidence level. The best fit parameters 
from the two dimensional likelihood function are fi^ = 0.09 and a* = 6.60. 

The one-dimensional likelihoods of individual cosmological parameters are shown in 



Fig. 5.4 The best fit values are a* = 2.3 and fi^ = 0.06. The la intervals are 0.3 < a < 8.9 
and 0.04 < < 0.14. 

At the moment lookback time versus redshift and Hubble parameter versus redshift 
measurements can not provide strong constraints on cosmological parameters. Both data 
sets are expected to improve significantly in the near future [8]j|- The updated data sets 
in combination with other measurements might prove to be very useful in constraining 
cosmological parameters. 
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Figure 5.3: 1, 2, and 3a confidence level contours on (j)CDM model parameters lookback 
time versus redshift data. Cross denotes the best-fit parameters at = 0.09 and a* = 0.60 
with = 20.54 for 30 degrees of freedom. 
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Figure 5.4: One dimensional likelihood functions of individual cosmological parameters 
from lookback time versus redshift data. The maximum likelihood values are normalized to 
one. 
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Chapter 6 



Constraints on ^CDM from SNela 
and GRB Data 

This chapter is based on Ref . [90l [91] . 

6.1 Supernova Type la as standard candles 

Supernovae (SNe) are highly luminous objects. At the peak of their activity SNe can be 
as bright as an average galaxy. They reach maximum brightness shortly after the stellar 
explosion and then slowly fade away in several weeks. The shock wave from explosion 
pushes a spherical shell of gas and dust into the interstellar medium and leaves behind a 
visible SNe remnant. SNe are classified based on the absorption lines of chemical elements 
in their spectra. SNe that do not have hydrogen atom Balmer lines in their spectra are 
classified as Type I, while those that have Balmer lines belong to type II. SN can be farther 
divided into subtypes. SNe type la (SNela) for example lack hydrogen lines but have a 
singly ionized silicon (Si II) line at 615 nm. SNela are the most useful for cosmological 
parameter estimation. First of all, they are astonishingly bright and can be seen at large 
distances. Secondly, they are very homogeneous in spectra, light-curve shapes, and peak 
absolute magnitude and can be easily calibrated and used as standard candles. 

SNela are believed to be carbon-oxygen white dwarfs in a binary system. The white 
dwarf accretes matter from its companion until it reaches the Chandrasekhar mass limit 
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of 1.4 solar masses. As the density and temperature inside the star increases, the mass 
in the core undergoes nuclear fusion, which results in a SNela explosion. The masses of 
white dwarfs at the time of explosion are very similar and this could explain why SNela 
have similar absolute magnitudes. The inhomogeneities in absolute magnitude of SNela are 
correlated with other observables. For example, brighter SNela tend to fade slower. This 
relationship can be used to further reduce the dispersion in absolute magnitude of SNela 
and make them a highly accurate standardizable candle. 

The measurement of absolute magnitudes of high-redshift SNela gave the first direct 
evidence of the accelerated expansion of the Universe (see Refs. |l6l[l7j). They were followed 
by more high-quality observations of SNela (see, e.g., Refs. [THH2ni ESHSII)- We now have 
more than 400 well calibrated, high-redshift SNela, going up to redshift of 1.7, that can be 
used to determine cosmological parameters. 

Ref. [3] used results of several new and old observations to compile a new SNela data set. 
They used a single, consistent analysis on all subsamples and implemented a new procedure 
to reject outliers. This compiled "Union" data set includes 307 SNela in the redshift range 
of z = 0.015 to 1.551. This data is consistent with the spatially-fiat ACDM model and gives 
best-fit value = 0.731+|J:[!^^(stat)^[!:J^^^(sys) in this model. The union data by themselves 
give more than 5a evidence for the accelerated expansion of the Universe. It has been also 
used to constrain parameters of other dark energy and modified gravity models (see, e.g.. 



6.2 Constraints on (/)CDM from Supernova Type la 



In astronomy, for historical reasons, flux (power per unit area) is usually given as apparent 



data 



magnitude 




(6.1) 
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where m{z) is the apparent magnitude of an object at redshift z, is its flux, and 0o is 
some arbitrarily chosen reference flux. Absolute magnitude M measures the objects intrinsic 
brightness and can be computed from apparent magnitude if the distance to the object is 
known. In cosmology fluxes are normalized so that the absolute magnitude of a SNela is 
equal to the apparent magnitude it would have if it were at a distance of 100 Kpc. Apparent 
magnitude then can be expressed through luminosity distance as 



m(z) = M + 25 + 51ogK(z)). (6.2) 
Redshifts, distance moduli fi = m — M, and measurement errors for SNela in the Union 



data set are given in Appendix D We can compute luminosity distance at different redshifts 
in the 0CDM model and compare theoretical predictions with the measurements. The 
luminosity distance depends on the parameter a, the nonrelativistic matter energy density 
Qrn, and the Hubble constant Hq. Hq can be determined from the data by minimizing the 
function, or it can be marginalized over with a prior probability distribution function 
PiHo). 



Constraints from SNela data on the 0CDM model are shown in Fig. 6.1 SNela data 
constrains the nonrelativistic matter density parameter to be less than 0.4 at about 3cr 
confidence level. The a parameter is constrained to be less than 4.5 at about 3a confidence 
level. 

The one dimensional likelihood functions for individual cosmological parameters are 



shown in Fig. 6.2 The best fit values are a* = and = 0.25. The la intervals are 



< a< 1.1 and < < 0.29. 

6.3 GRB as standard candles 

One way of improving our understanding of how dark energy behaves is to study the evolu- 
tion of the Universe at redshifts higher than those probed by SNela. This requires standard 
candles that are visible at greater distances. Gamma-ray bursts (GRBs) could in principle 
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Figure 6.1: 1, 2, and 3a confidence level contours on (f) CD M model parameters from SNela 
absolute magnitude versus redshift data. Cross denotes the best-fit value at fij^ = 0.29 and 
a* = with = 320 for 305 degrees of freedom. 

serve as such high redshift standardizable candles. They are the most luminous events in 
the Universe today and can be seen to beyond z = 8. If it is definitely established that 
GRBs are standardizable candles, their visibility at high redshift should prove to be very 
useful in discriminating between ACDM and time-varying dark energy models. 

With the intention of getting cosmological constraints from GRB observations a number 
of GRB calibrations have been used so far [HS]. One that gives least scatter and therefore 
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Figure 6.2: One dimensional likelihood functions of individual cosmological parameters 
from SNela absolute magnitude versus redshift data. The maximum likelihood values are 
normalized to one. 

most information is 

a relation that connects the total burst energy of the GRB {E^) to the peak energy of the 
GRB spectrum (-Epeak) [96j. Regrettably, we do not yet have a model-independent way of 
computing the coefficients Ai and Bi. A better understanding of physical processes that 
result in the burst, or observations of nearby GRBs (to which distances can be measured 
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independently), could in principle help us to calibrate the E^-Epeak relation without any prior 
assumptions. Until recently the only way to extract cosmological information from GRBs 
was to recalibrate them for every dark energy model considered (at each set of parameter 
values). This is time consuming and also results in large statistical uncertainties and hence 
GRB cosmological constraints that are poor. 

Recently, two methods of calibrating GRBs in cosmology-independent manners have been 
proposed and used to constrain some dark energy models. One method externally calibrates 
GRBs, by using SNela measurements pTlllOU] . while the other uses GRB measurements 
alone |101] . While externally calibrating GRBs results in tighter cosmological constraints, 
internally calibrated GRB data may be straightforwardly combined with other data when 
deriving cosmological constraints. The resulting cosmological constraints are still loose, 
but in the future when more high precision GRB observations become available this could 
provide a strong test of dark energy. 

Reference |lUlj recently used data of 69 GRBs [HS] to construct a model-independent 
distance measure that can be used to constrain cosmological models. When this method is 
used to constrain ACDM the GRB data favor lower values of both cosmological constant 
energy density (^^a) and nonrelativistic matter energy density (i^m) than do the SNela data. 
The GRB data by themself are unable to strongly constrain cosmological parameters, for 
example in spatially-fiat ACDM the GRB data require fim = 0.25^q at la confidence 

For early discussions of the use of GRBs as a cosmology probe see, e.g., Refs. [96l I102f 
I105j . More recent studies may be traced back through Refs. |106m09j . For a review of GRB 
physics see, e.g., Ref. |110] . 

6.4 Constraints from GRB data 



Besides the Epeak-E^ relation, Eq. (6.3), four other calibrations for GRBs are used in 



Ref. |101j . Reference |101] placed each of the 69 GRBs in the redshift range z = 0.17 
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to z = 6.6 at a luminosity distance that minimized a combined thai took weighted ac- 
count of all five calibration relations. From this |lUlj computed the cosmology-independent 
distance measure 



where 



Ho 1 



(6.4) 



(6.5) 



hcz{i + zy/^ 

and dilz) is the luminosity distance at redshift z, and c is the speed of light. The ratio in 



Eq. (6.4) does not depend on the Hubble constant and does not require information about 



the absolute calibration of GRBs (which are unknown). 

Reference |1U1] computed the cosmology-independent distance measure fp in six redshift 
bins fp{zi), i = 1, 2, ... 6. The values of fp{zi) are shown in Appendix [E] and the normalized 
covariance matrix is 



S 



(6.6) 



/ 1.0000 0.7056 0.7965 0.6928 0.5941 0.5169 \ 

0.7056 1.0000 0.5653 0.6449 0.4601 0.4376 

0.7965 0.5653 1.0000 0.5521 0.5526 0.4153 

0.6928 0.6449 0.5521 1.0000 0.4271 0.4242 

0.5941 0.4601 0.5526 0.4271 1.0000 0.2999 

\ 0.5169 0.4376 0.4153 0.4242 0.2999 1.0000 / 

For currently viable cosmological models, these rp{zi) are almost completely independent of 
the cosmological model and so provide a useful summary of current GRB data |101] . This 
information can be used to constrain any dark energy model and the resulting GRB data 
constraints can be straightforwardly combined with other constraints. 

We compute the difference between the theoretical prediction and the measured value 
at each of the 6 redshifts, 



r^rdata 
P 



(;,.)_ ^theory 



and the function 



x'(Q^,p) = A(^O^i(^-')ij^jA(^j) 



(6.7) 



(6.8) 



42 



where 5*^ is the normahzed covariance matrix given by Eq. (6.6) and summation over re- 
peated indexes is assumed. Here a\ is if ^{z^) > and cr^ if A{zi) < 0, see Appendix [Ej 



The constraints on (j)CDM from GRB data are shown in Fig. |6.3[ Current GRB data can 
not constrain cosmological parameters well and a wide range of values are at 3a confidence 
level. The best-fit parameter values are a* = 10.2 and = 0.0 with = 1-39 for 4 
degrees of freedom. 




Figure 6.3: 1, 2, and 3a confidence level contours on (pCDM model parameters from GRB 
data. Cross denotes the best-fit value at VL*^ = 0.29 and a* = 0. 
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The one dimensional likelihood functions for individual parameters are shown in Fig. 6.4 
The best fit parameter values are a* = 2.3 and fij^ = 0. One standard deviation intervals 
are < a < 12.3 and < < 0.17. 

1i ^ ^ 1 1 




a 




Figure 6.4: One dimensional likelihood functions of individual cosmological parameters 
from GRB data. The maximum likelihood values are normalized to one. 

GRB data by themselves do not provide tight constraints on cosmological parameters. 
Moreover, while not greatly significant, current GRB data favor cosmological parameter 
values that are at odds with what other data favor. When used in combination with some 
of the highest quality current data (e.g., SNela and BAO peak measurements) current GRB 
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data only slightly change the results. This is mainly because in the absence of an absolute 
calibration of GRBs they, as standard candles, have big measurement uncertainties. This 
is however quite likely to change as more and better-quality GRB measurements become 
available and with improvements in how to calibrate GRBs. GRBs could potentially provide 
a very strong test of the time variation of dark energy as they can be observed up to redshifts 
beyond 8, at distances where other standard candles cannot be detected. 
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Chapter 7 

Constraints from Cluster Gas Mass 
Fraction Data 

This chapter is based on Ref. 

7.1 Using galaxy clusters as cosmological probes 

Galaxy clusters are the biggest gravitationally bound objects in the Universe. Old, relaxed, 
rich ones should provide a fair sample of the matter content of the Universe. The ratio 
of baryonic mass to total nonrelativistic mass in the clusters should be close to the ratio 
of cosmological parameters i^b/^m, the baryonic fraction /. More than 80% of clusters' 
baryonic mass is in hot X-ray emitting intergalactic gas. Estimates of gas mass fraction in 
clusters can be combined with the priors on Q^, from WMAP or BBN to get tight constrains 
on nonrelativistic matter density. During cosmic evolution the ratio Q^/Qra should remain 
the same, but clusters are observed at different redshifts and since the reconstructed /gas 
depends on the assumed distance to the cluster, this data can be used to constrain cosmic 
evolution and hence the dark energy. 

To get a good estimate of /gas, the gas mass fraction, large, relaxed clusters in thermal 
equilibrium should be used. Clusters with high core temperature {kT > 5 keV) are better 
for the analysis, because in computer simulations and observations they have less systematic 
scatter in gas mass fraction compared to "cold" clusters. 



46 



The gas mass fraction is given by |6] 



KA^bjz) Ufjz) 

1 + S{z) V dA{z) 



) 



1.5 



/gas(^) 



(7.1) 



Here is a calibration constant that accounts for effects, such as accuracy of instrument 



cahbration and X-ray modelhng. A conservative 10% Gaussian uncertainty K = 1.0 ± 0.1 
is used in Ref. f6j. The factor A accounts for the change in the angle subtended by the 
cluster between the reference cosmology and the cosmology of interest and is always very 
close to 1. 7 models non-thermal pressure due to subsonic motion in the gas, cosmic rays, 
magnetic fields, etc. In the analysis of Ref. [6] uniform prior with 1.1 < 7 < 1.2 is assumed. 
b is the "bias" factor, the ratio by which the baryon fraction in the cluster differs from the 
same mean ratio in the Universe. The "bias" factor is modeled as 6 = 6o(l + Q^b^) [6] and 
cosmological simulations suggest the uniform priors 0.65 < 60 < 1-0 and —0.1 < ab < 0.1, 
[6]. The parameter s = Sq{1 + agz) models the baryon gas mass fraction in stars. We use 
the uniform prior with —0.2 < Os < 0.2 and the Gaussian prior with sq = 0.13 ± 0.01, [6]. 
(i^^ is the angular diameter distance computed in a reference, spatially-fiat ACDM model 
with Q\ = 0.7, and is the angular diameter distance computed in the 0CDM model. 

7.2 Many dimensional integration 

We use the modified technique discussed below to perform some of the integrals to account 
for these nuisance parameters, which have usually veen done numerically, as discussed next. 

Let's say we want to integrate a function f{x) on the interval (a, 6). If this cannot be 
done analytically, different numerical integration algorithms may be used to compute the 
integral to a given precision. 

If the problem is one dimensional and the integrand is reasonably well-behaved, (a, b) can 
be subdivided into segments of equal length; in each segment /(x) can be interpolated by 
a polynomial; and the integral can be approximated by a sum |112] . The simplest method 
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of this sort uses the trapezoidal Newton-Cotes formula 



The error on the integral computed using Eq. (7.2) can be approximated by 



^i~^sup|/'(x)| (7.3) 

where prime denotes the derivative. The error could be smaller than this if the function is 
sufficiently smooth. 

This method can be applied to multi-dimensional integrals also, by transforming them 
into repeated one-dimensional integrals. The number of function evaluations required to get 
the same precision, however, will grow exponentially in this case and so will the required 
computational time. To sample a unit interval in one dimensional space with a distance 
between adjacent points of 0.01, only hundred evenly-spaced points are required. To do the 
same thing in 10 dimensional space will require lO^'' points. 

To cut down on computational time Monte-Carlo integration techniques can be em- 
ployed. Monte-Carlo algorithms randomly pick points in the multi-dimensional space. 
By the law of large numbers, the error on an integral computed using Monte-Carlo method 
scales as 1/ independent of dimensionality. 



The gas mass fraction in Eq. (7.1) depends on eight "nuisance" parameters. These are 
7, 6o, «b, So, Os, ^b, and Hq. To perform this eight-dimensional integration we use a 
method that in this particular case works faster than Monte-Carlo sampling. 



Five out of eight nuisance parameters enter Eq. (7.1) as a redshift independent combina- 
tion r = Kjbo{nbh'^)/h^^^. We can think of each nuisance parameter as a random number 
with a given probability distribution function (PDF). Then F will also be a random variable 
with some PDF. 

We use the following statistics results |113] . If two random variables a and b are inde- 
pendent with PDFs Pa{x) and Ph{x) then variables c = ah, d = a/b and / = F{a) are also 
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random with PDFs 

1 




Pc= P^{x')Pk{x")6{x'x" - x)dx'dx" = / —P^{x')P^{x/x')dx', (7.4) 



\x 



~ J J Pa.{x)Ph{x")5{x / x" — x)dx dx" = j \x\Pa.{xx)P\,{x)dx\ (7.5) 



^^'^"^^ P.{F'\x)). (7.6) 



dx 



We first numerically compute the PDF for P using Eqs. (7.4 7.6) and replace the the five-fold 
integration of the likelihood function by a one- dimensional integral. This reduces computa- 
tional time significantly. 

7.3 Constraints from /gas data 

We use the measurements of gas mass fraction of 42 hot, X-ray luminous, relaxed galaxy 
clusters in the redshift range oi z = 0.01 to 1.1 [6j. The cluster gas mass fraction data is 
presented in Appendix [F} We integrate the likelihood over the Hubble constant with the 
Gaussian WMAP prior h = 0.73 ± 0.03, and baryonic mass energy density with a Gaussian 
prior fib = (0.0223 ± O.OOOS)/^-^. 



The constraints on 0CDM parameters are shown in Fig. |7.1[ Brown, light blue and deep 
blue areas correspond to 1, 2, and 3a confidence levels respectively. The black cross denotes 
best-fit parameters fijjj = 0.27 and a* = 0.0, with = 43.5 for 40 degrees of freedom. 
The nonrelativistic matter density is constrained to be in the range of 0.15 < < 0.45 
at about 3cr confidence level and a is constrained to be less than 5 at about 3a confidence 
level. 



One dimensional likelihood functions for individual parameters are shown in Fig. 7.2 
The best-fit values from the one dimensional likelihood functions are a* = 0.0 and = 
0.27. la confidence level intervals from the one dimensional likelihood functions are < 
a < 1.1 and 0.22 < Qn. < 0.33. 
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Figure 7.2: One dimensional likelihood functions of individual cosmological parameters 
from galaxy cluster gas mass fraction data. The maximum likelihood values are normalized 
to one. 
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Chapter 8 

Constraints on ^CDM from Baryon 
Acoustic Peak Measurements 

This chapter is based on Refs. [901 1114] . 

8.1 Transverse baryon acoustic peak measurements 

Before recombination baryons and photons are tightly coupled and gravity and pressure 
gradients induce sub-acoustic-Hubble-radius oscillations in the baryon-photon fluid |115[ 
I116j . These transmute into the acoustic peaks observed now in the CMB anisotropy angular 
power spectrum, which provide very useful information on various cosmological parameters. 
The baryonic matter gravitationally interacts with the dark matter and so the matter power 
spectrum should also exhibit these "baryon acoustic" wiggles. Because the baryonic matter 
is a small fraction of the total matter the amplitudes of the BAO wiggles are small. The BAO 
peak length scale is set by the sound horizon at decouphng, ~ 10^ Mpc, and so detecting 
the BAO peak in a real space correlation function requires observationally sampling a large 
volume. The BAO peak in the galaxy correlation function has recently been detected by 
using SDSS data [H 1117] and by using 2dFGRS data |118j . For more recent discussions of 
the observational situation see Refs. [HI I119[ I12U] . 

The sound horizon at decoupling can be computed from relatively well-measured quanti- 
ties by using relatively well-established physics. Consequently it is a standard ruler and can 
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be used to trace the universe's expansion dynamics |12 11 - 1124] . A measurement of the BAO 
peak length scale at redshift z fixes a combination of the angular diameter distance and 
Hubble parameter at that redshift. More precisely, what is determined [7] is the distance 

Dy{z)={{\^zfd\{z)zlH[z)\^l\ (8.1) 

where li{z) is the Hubble parameter and dA the angular diameter distance. Dy{z) depends 
on the cosmological parameters of the model, including those which describe dark energy, so 
we can constrain these parameters by comparing the predicted Dyi^z) to the measurements. 

We examine the constraints on 0CDM from two measurements of the BAO peak. The 
first is from the BAO peak measured at z = 0.35 in the correlation function of luminous 
red galaxies in the SDSS This measurement results in A{0.35) = 0.469 ± 0.017 (one 
standard deviation error), where the dimensionless and ifo-independent function 

A(z) = Dv^(z)V^5S ^82) 



and Dv{z) is the distance measure defined in Eq. (8.1). The measured value of A(0.35) 
does not depend on the dark energy model and only weakly depends on the baryonic energy 
density. The measurement also has a weak dependence on parameters like the spectral index 
of primordial scalar energy density perturbations (the assumed value is n = 0.98) and the 
sum of the neutrino masses, but this is not strong enough to have significant effect on the 
final result. To constrain cosmological model parameters in this case we perform a standard 
analysis. 

The second BAO peak measurement we use is from the correlation function of galaxy 
samples drawn from the SDSS and 2dFGRS at two different redshifts, z = 0.2 and z = 0.35, 
as determined by Ref. [Sj. This analysis includes the SDSS luminous red galaxies, so the 
the two BAO peak measurements are not statistically independent. 

This measurement gives the correlated values Vg/ Dv{0.2) = 0.1980±0.0058 and rs/Dv{0.35) 
0.1094 ± 0.0033 (one standard deviation errors), where is the comoving sound horizon at 
recombination. These two measurements are correlated, with the inverse of the correlation 
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matrix given b 



V 



-1 



35059 -24031 
-24031 108300 



To compute we first compute the angular diameter distance to the surface of last scatter- 
ing, ^^(1089). We then use the WMAP measurement of the apparent acoustic horizon angle 
in the CMB anisotropy data |125j to determine the sound horizon = [(1 + 2;)c?a(-2)] |2=io89 x 
0.0104 (where we ignore the WMAP measurement uncertainty and assume that is known 
perfectly). The use of the WMAP prior on the apparent acoustic horizon angle results in 
very tight constraints on the spatial curvature. When this measurement is not used, these 
measurements alone can not tightly constrain the dark energy parameters. 

To constrain cosmological parameters in this case we follow Ref. P] and first compute 

/ r,/Dy(0.2,f]^,«)- 0.1980 \ 
l^^m, «j - r,/Dy(0.35, fin,, «) - 0.1094 ) ' ^^-^^ 

where for definiteness we consider the 0CDM model. We then compute the function 

x\nn.a)=X-^V-'X. (8.4) 

and the likelihood function 

L(fi^, a) oc exp(-x'(fim, a)/2). (8.5) 

8.2 Constraints from transverse BAO peak measurem- 
nts 



The constraints on (/)CDM model from BAO data presented in Ref. [7j are shown in Fig. 8.1 



Constrains from BAO data presented in Ref. [8] are shown in Fig. 8.2 The one dimensional 



individual likelihoods are shown in Figs. 8.3 and 8.4 respectively. The best fit values of 



individual parameters from Fig. 8.4 are a* = and = 0.24, with la intervals < a < 4.5 
and 0.24 < fi„, < 0.26. 

These measurements were made at only two redshifts and by themselves do not provide 
a robust and reliable test of dark energy models, but in combination with other data they 
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do provide useful constraints on cosmological parameters. In addition, a number of surveys 
are planned in next few years that will measure the BAO scale accurately and at a variety 
of redshifts up to z = 1.2. This upcoming BAO data, especially when combined with other 
data, will prove very useful in tightly constraining dark energy parameters [126j. 




Figure 8.1: 1, 2, and 3a confidence level contours on (pCDM model parameters from BAO 
peak measurement of Ref. 
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0.35 



Figure 8.2: 1, 2, and 3a confidence level contours on (pCDM model parameters from BAO 
peak measurement of Ref. f^. Cross denotes the best-fit value at = 0.25 and a* = 0. 

8.3 Radial BAO peak measurements 

Recently it was argued tliat tlie above measurements of tlie BAO scale were essentially 
measurements orthogonal to the line of sight and so statistically independent from a line 
of sight measurement of the BAO scale, even if the same galaxy catalog is used for both 
measurements. Reference jl27j used the SDSS data to compute the line of sight or radial 
BAO scale in redshift space for two ranges of redshift and showed the resulting constraints 
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Figure 8.3: One dimensional likelihood functions of individual cosmological parameters 
from BAO peak measurement of Ref. The maximum likelihood values are normalized to 
one. 



on the spatially- flat XCDM parameterization [128J. These constraints from the radial BAO 

scale data are quite similar to the constraints derived from earlier "transverse" BAO scale 
measurements. The data are consistent with spatially-fiat ACDM. However, these current 
radial BAO measurements (like current non-radial BAO measurements) can not tightly 
constrain time-varying dark energy by themselves (although the situation is anticipated to 
improve in the next few years), as discussed next. 
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Figure 8.4: One dimensional likelihood functions of individual cosmological parameters 
from BAO peak measurement of Ref. 18J. The maximum likelihood values are normalized to 
one. 

In a spherically symmetric Universe the two-point correlation function is a function of two 
variables, ,^ = ^(cr, tt), where a is the separation along the line of sight and vr is the separation 
on the sky. It can also be expressed as a function of absolute separation r = ^/a^ + n'^ and 
the cosine of the angle between the line of sight and the direction of separation, /i = vr/r. 
The correlation function can then be decomposed into multipole moments, 

a(r) = ^^ ar,fi)Pi{fi)dfi, (8.6) 
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where Pi is the /* order Legendre polynomial. Multipole moments of different orders can 
be related to each other if one has a complete theory of linear and nonlinear evolution. 
Although high multipoles that describe the "shape" of baryon acoustic oscillation imprints 
on the matter distribution are very difficult to measure in practice, theoretically they are 
independent of the monopole and could provide additional structure formation tests. 

Initial work considered only the averaged over direction monopole part of the correlation 
function, 

Ur) = l ar,f^)dfi, (8.7) 

and found a BAO peak signal at a comoving distance of r ^ 110/i~^Mpc. This measurement 
was however mostly transverse to the line of sight direction vr; the weight of separation along 
the line of sight contributes less then 1%. Consequently, it is fair to assume that the radial 
baryon acoustic peak scale measurement in the line of sight direction from ^ (cr) is statistically 
independent from that measured from ^(r) ~ ^(^r), even if the same galaxy sample is used 
for both measurements. 

Reference. |127] used SDSS data to measure the radial baryon acoustic scale in two 
redshift ranges 2 ~ 0.15 - 0.30 with radial BAO peak scale Az = 0.0407 ± 0.0014 and z ~ 
0.40 — 0.47 with Az = 0.0442 ± 0.0016 (both one standard deviation errors). Theoretically 
the radial BAO peak scale is given by 

Az = H{z)rsizd)/c (8.8) 

where H{z) is the Hubble parameter at redshift z, rs{zd) is the sound horizon size at the drag 
redshift z^, at which baryons were released from photons, and c is the speed of light. H{z) 
can be easily computed in a given cosmological model and depends on model parameters 
such as the non-relativistic matter density and the time dependence of dark energy. 

Ts can be computed with two different methods. One is to use the ratio Ig between the 
distance to the last-scattering surface and measured by CMB anisotropy experiments and 
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compute the sound horizon at photon decouphng from 



= + (8.9) 



''S 



Here z=k is the redshift at photon decouphng and d-A is the angular diameter distance. Aher- 
natively, one can use priors on the fractional energy density parameters of baryonic matter, 
fib; nonrelativistic matter, Q^, and relativistic matter, Qj., from, e.g., CMB anisotropy 
measurements, and compute the sound horizon at the drag redshift from 

r,{z^) = ^ / =. (8.10) 

Ho^/m^J^ ^{a + 1.69fir/fim)(l + a0.75fib/fir) 

Both options have similar drawbacks. One has to assume priors on "nuisance" parameters 
like Is or various energy densities. CMB anisotropy measurements themselves have measure- 
ment errors that must be accounted for, otherwise the errors on the estimates of dark energy 
model parameters of interest will be underestimated. Also, the best fit values for nuisance 
parameters given by CMB anisotropy data are in general different for every cosmological 
model and also depend on model parameter values. To be fully consistent when using priors 
one would have to reanalyze CMB experiments for each cosmological model (and model 
parameter value) instead of using a single set of values for Ig, fib, fim, and fir. 

At present, however, the BAO scale is measured only in two redshift ranges and does 
not provide very tight parameter constraints compared to other observational tests. Hence, 
as long as we are interested in preliminary constraints on dark energy from BAO scale 
measurements we may use the simplified approach of Ref . |127] , keeping in mind that when 
more and better quality BAO scale measurements become available a more complete, careful, 
and time-consuming analysis will be warranted. 

8.4 Constraints from radial BAO peak measurements 



The confidence level contours for the spatially-flat </)CDM model are shown in Fig. 8.5 



and 8.6 Corresponding one dimensional likelihood functions for individual cosmological 
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parameters are shown in Fig. 8/7 and 8^ Here, the radial BAO measurements constrain 
to be between 0.15 and 0.4 at about 3a, but the a parameter is not constrained well 
and large values of a (relatively rapidly evolving dark energy) are not ruled out, although 
the likelihood peaks at a = 0. These results are similar to the ones derived earlier using the 
non-radial BAO peak scale measurements. 




Figure 8.5: 1, 2, and 3a confidence level contours on (pCDM model parameters from radial 



BAO peak measurements using the technique of Eq. (8.9). The best fit values are a* = 
and = 0.23. 
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Figure 8.7: One dimensional likelihood functions of individual cosmological parameters 



from radial BAO peak measurements and using the technique of Eq. (8.9). The maximum 
likelihood values are normalized to one. 
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Figure 8.8: One dimensional likelihood functions of individual cosmological parameters 
from radial BAO peak measurements and using the technique of Eq. (8.10). The maximum 
likelihood values are normalized to one. 
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Chapter 9 
Conclusion 



9.1 Joint constraints 

Chapters [5] to |8] present constraints on (pCDM model parameters using different available 
cosmological data sets. The lookback time versus redshift and the GRB data cannot at the 
moment constrain dark energy strongly. Likelihoods resulting from these two data sets are 
flat in the direction of a and give preference to cosmological models with low nonrelativistic 
mass energy density < 0.2, which is at odds with other observations, but this is not very 
significant. 

In both cases this happens because the measurements are available only at a small 
number of redshifts (32 for lookback time data and 6 for GRB) and the measurement 
uncertainties are large. This should change when more and better quality data become 
available. In the future these two methods could prove to be very useful in constraining 
dark energy. 

Constraints from SNela absolute magnitude versus redshift data, galaxy cluster gas 
mass fraction versus redshift data, and BAO peak measurements are considerably tighter. 
BAO peak measurements are not very effective in constraining the a parameter at the 
moment, but they give very good constraints on nonrelativistic matter energy density. The 
power of the constraints seems surprising given that the BAO peak measurements have been 
made at only two redshifts so far. SNela and cluster gas mass fraction data are not that 
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effective in constraining but provide upper limits on a. Current data favors tlie time- 
independent cosmological constant, but tliere is a large part of time-evolving dark energy 
((/>CDM) parameter space that is still not ruled out. 

Since none of the available data sets alone can simultaneously tightly constrain both cos- 
mological parameters a and Q^n, we derive tighter constraints by performing joint analyses 
of observational data. We get the best results when using SNela and BAO peak measure- 
ments; adding other data to the mix does not significantly improve the constraints. Since 
these two are independent we define the joint likelihood as the product of the individual 
likelihoods 

Aot = '^SN'^BAO- (9-1) 

The constraints derived from this joint likelihood function are shown in Fig. |9.1[ The joint 
analyses results in much tighter constraints. The nonrelativistic matter density is in the 
range of 0.18 < < 0.32 and a is constrained to be less than 1.5 at about 3a confidence. 
The best-fit values a* = 0.27 and = 0.24 correspond to slowly-varying dark energy but 
the time-independent cosmological constant also fits the data very well. 



The likelihood functions of individual cosmological parameters are shown in Fig. 9.2 
The best-fit values from one dimensional likelihood functions are a* = 0.27 and fl^ = 0.24. 
The la intervals are 0.02 <a< 0.54 and 0.22 < < 0.26. 

9.2 Future prospects of detecting time-evolving dark 
energy 

Joint analysis of currently available cosmological data already results in strong constraints 
on the (/)CDM model, fl^i is estimated to be in between 0.2 and 0.3 and a is constrained to 
be less than 1.5. 

The quality and quantity of available cosmological data sets is expected to increase 
significantly in the next decade. This will result in much tighter constraints on 0CDM 
model parameters. Some of the uncertainty in a and comes from the fact that we do 
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0.15 



0.2 



0.25 



0.3 



0.35 



Figure 9.1: 1, 2, and 3a confidence level contours on (pCDM model parameters from joint 
analysis of SNela Union data and transverse BAO peak measurements. The best-fit param- 
eters are a* = 0.3 and fl^ = 0.24 with = 321 for 307 degrees of freedom. 

not know precise values of nuisance parameters, such as Hq and fib- Ongoing surveys (such 
as HST) and new missions (such as PLANCK satellite) will measure these with greater 
accuracy, which will translate into better constraints on a and Q^. 

From the joint analysis of near future SNela, GRB, galaxy cluster gas mass fraction and 
BAO measurements we should be able to either detect the time-dependence of dark energies 
energy density at a high confidence level or to constrain it to small, physically uninteresting 
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Figure 9.2: One dimensional likelihood functions of individual cosmological parameters 
from the joint analysis of SNela Union data and transverse BAO peak measurements. The 
maximum likelihood values are normalized to one. 



values. 
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Appendix A 

CH — h function for a numerical 
solution of the scalar field equations. 



#include<iostream> 
#include<fstream> 
#include <cmath> 

using namespace std; 



double fy( double a, double y, double dy , double kmp, double alpha); 
double fa (double a, double y, double dy , double kmp, double alpha); 



// 

// This function takes as an input the value of alpha, nonrelativist ic 
// matter density Omegajn, and redshift z. 

// It solves ODEs using 4th order Runge— Kutta method and returns 
// comoving distance. 

// 



double dist ance ( double alpha, double OmegaM, double z) 
{ 

const double tO = 0.01; 
const double dt = 0.001; 

double t = to ; 

double A= 2 . /3 . * alpha *( alpha + 2.0); 

double kmp = 8 . 0/3 . * ( alpha + 4.0)/(alpha + 2.0)*pow(A , alpha/2.0); 



double afin , dafin ; 
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double sum = 0.0; 

// Initial conditions on scale factor — a, and scalar field — y. 
double a = pow(tO , 2.0/3.0); 

double y = sqrt (A) *pow ( tO , 2.0/(alpha + 2.0)); 
double dy = 

sqrt (8.0/3.0* alplia/(alpha + 2.0))/pow(t0 , alpha /( alpha + 2.0)); 
double OM = 

4.0/9.0/a/a/a/(4.0/9.0/a/a/a + (dy*dy + kmp/pow(y , alpha ) ) / 1 2 . ) ; 

double kll, kl2 , kl3 , k21 , k22 , k23 , k31 , k32 , k33 , k41 , k42 , k43 ; 

while (OM > omegam) 
{ 



kll = 


fa (a, y, dy , kmp, alpha )*dt; 










kl2 = 


dy* dt ; 










kl3 = 


fy(a, y, dy , kmp, alpha)*dt; 










k21 = 


fa(a + 0.5*kll, y + 0.5*kl2, dy + 0.5*kl3, 


kmp. 


alpha ) 


* dt 


k22 = 


(dy + 0.5*kl2)*dt ; 










k23 = 


fy(a + 0.5*kll, y + 0.5*kl2, dy + 0, 


.5*kl3 , 


kmp. 


alpha ) 


* dt 


k31 = 


fa(a + 0.5*k21, y + 0.5*k22, dy + 0, 


.5*k23 , 


kmp. 


alpha ) 


* dt 


k32 = 


(dy + 0.5*k22)*dt ; 










k33 = 


fy(a + 0.5*k21, y + 0.5*k22, dy + 0, 


.5*k23 , 


kmp. 


alpha ) 


* dt 


k41 = 


fa (a + k31 , y + k32 , dy + k33 , kmp. 


alpha ) 


* dt ; 






k42 = 


(dy + k32)*dt ; 










k43 = 


fy(a + k31 , y + k32 , dy + k33 , kmp. 


alpha ) 


* dt ; 







a += kll/6.0 + k21/3.0 + k31/3.0 + k41/6.0; 
y += kl2/6.0 + k22/3.0 + k32/3.0 + k42/6.0; 
dy += kl3/6.0 + k23/3.0 + k33/3.0 + k43/6.0; 



t = t + dt ; 

0M= 4.0/9.0/a/a/a/(4.0/9.0/a/a/a + (dy*dy + kmp/pow(y , alpha ) ) / 1 2 . ) ; 

} 

af in = a ; 

dafin = fa (a, y, dy , kmp, alpha); 
t = to ; 

a = pow(tO , 2.0/3.0); 

y = sqrt (A)*pow(tO , 2.0/(alpha + 2.0)); 
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dy = sqrt (8.0/3. 0*alpha/(alpha + 2.0))/pow(t0 , alpha/(alpha + 2.0)); 
OM = 4.0/9. 0/a/a/a/(4. 0/9. 0/a/a/a + (dy*dy + kmp/pow(y , alpha ) ) / 1 2 . ) 

while (OM > omegam) 
{ 

kll = fa(a, y, dy, kmp, alpha )*dt; 

kl2 = dy*dt ; 

kl3 = fy(a, y, dy , kmp, alpha)*dt; 

k21 = fa(a + 0.5*kll, y + 0.5*kl2, dy + 0.5*kl3, kmp, alpha)>Kdt; 

k22 = (dy + 0.5*kl2)*dt ; 

k23 = fy(a + 0.5*kll, y + 0.5*kl2, dy + 0.5*kl3, kmp, alpha)*dt; 

k31 = fa(a + 0.5*k21, y + 0.5*k22, dy + 0.5*k23, kmp, alpha)*dt; 

k32 = (dy + 0.5*k22)*dt ; 

k33 = fy(a + 0.5*k21, y + 0.5*k22, dy + 0.5*k23, kmp, alpha)*dt; 

k41 = fa(a + k31 , y + k32 , dy + k33 , kmp, alpha)*dt; 

k42 = (dy + k32)*dt ; 

k43 = fy (a + k31 , y + k32 , dy + k33 , kmp, alpha )*dt; 

a += kll/6.0 + k21/3.0 + k31/3.0 + k41/6.0; 
y += kl2/6.0 + k22/3.0 + k32/3.0 + k42/6.0; 
dy += kl3/6.0 + k23/3.0 + k33/3.0 + k43/6.0; 

t = t + dt ; 

OM = 4.0/9.0/a/a/a/(4.0/9.0/a/a/a + (dy*dy + kmp/pow(y , alpha))/12. 

if(a/afin > 1.0/(1.0 + z)) 
{ 

sum += 1 . / a ; 

} 

} 

sum *= dt ; 
sum *= dafin ; 



return sum; 

} 

double fy( double a, double y, double dy , double kmp, double alpha 
{ 

return ( -3.0* sqr t (4 . / 9 . / a/a/a + 

1 . 0/ 1 2 . 0* ( dy*dy + kmp/pow((y) , alpha)))* 
dy + kmp* alpha /2 . 0/pow (( y) , alpha + 1.0)); 

} 
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double fa(double a, double y, double dy , double kmp, double alpha) 
{ 

return (sqrt(4.0/9.0/(a) + (a)*(a)/12.0*(dy*dy + kmp/pow(y , alpha ) ) ) ) ; 

} 
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Appendix B 

Ages of passively evolving galaxies. 



Table B.l: Ages of passively evolving galaxies versus red- 
shift data. 



z 


Age (Gyr) 


0.117100 


10.2000 


0.117400 


10.0000 


0.222000 


9.00000 


0.231100 


9.00000 


0.355900 


7.60000 


0.452000 


6.80000 


0.575000 


7.00000 


0.644000 


6.00000 


0.676000 


6.00000 


0.833000 


6.00000 


0.836000 


5.80000 


0.922000 


5.50000 


1.17900 


4.60000 


1.22200 


3.50000 


1.22400 


4.30000 


1.22500 


3.50000 


1.22600 


3.50000 


1.34000 


3.40000 


1.38000 


3.50000 


1.38300 


3.50000 


1.39600 


3.60000 


1.43000 


3.20000 


1.45000 


3.20000 


1.48800 


3.00000 


Continued on next page 
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Table B.l — continued from previous page 



z 


Age (Gyr) 


1.49000 


3.60000 


1.49300 


3.20000 


1 f>i nnn 


9 snnnn 


1.55000 


3.00000 


1.57600 


2.50000 


1.64200 


3.00000 


1.72500 


2.60000 


1.84500 


2.50000 
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Appendix C 

Hubble parameter versus redshift 
data 



Table C.l: Hubble parameter versus redshift data with 
1(7 measurement uncertainties. 



z 


H{z) 




0.09 


69 


12 


0.17 


83 


8.3 


0.27 


70 


14 


0.4 


87 


17.4 


0.88 


117 


23.4 


1.3 


168 


13.4 


1.43 


177 


14.2 


1.53 


140 


14 


1.75 


202 


40.4 
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Appendix D 

SNela "union" data set 



Table D.l: SNela "union" data set. The redshift dis- 
tance modulii and Icr statistical measurement errors 
on the measurement of 



z 






0.050043 


36 


7730018543 





153148680782 


0.052926 


36 


8241423261 





148454214927 


0.02513 


35 


1678175291 





165420933733 


0.070086 


37 


567740156 





15403068043 


0.062668 


37 


5515804571 





145317136417 


0.087589 


38 


1915656352 





159495013964 


0.078577 


37 


5212710132 





147429338762 


0.017227 


34 


6546665523 





1880259565 


0.042233 


36 


3586532798 





152312134352 


0.045295 


36 


6566139461 





148758275509 


0.019599 


34 


5182384466 





179081303451 


0.100915 


38 


5088119689 





145279087383 


0.027342 


35 


2300574166 





171939536277 


0.074605 


37 


6663645268 





147931348064 


0.026489 


35 


5199700817 





182864910719 


0.049922 


36 


6999653546 





146530609404 


0.030604 


35 


5404531616 





161624590233 


0.016345641 


34 


0924538801 





144775220876 


0.0154363 


34 


0329851021 





150597548663 


0.030529 


35 


5818149937 





0920067053122 


0.024525 


34 


9577174302 





109085696057 


0.023953 


34 


935520479 





115625212493 


0.026038 


35 


3566237332 





106784897071 


0.048948 


36 


6391120177 





1766982111 


Continued on next page 
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Table D.l - continued from previous page 



z 


11 




0.024314 


35 


0794913195 


0.18643839963 


0.015166 


34 


2041114924 


0.225518080121 


0.03572 


36 


1050539158 


0.188094161827 


0.048818 


36 


3507343053 


0.166961986591 


0.0219800059146 


34 


8811235819 


0.196258768364 


0.0275 


35 


6735856723 


0.181507304221 


0.1244 


38 


9521645786 


0.168017789501 


0.036 


35 


8260342723 


0.173778570484 


0.01673 


34 


1592687583 


0.209352704395 


0.016321 


34 


1144805357 


0.209884301455 


0.021793 


34 


9048499767 


0.282608881266 


0.01645 


34 


1818298266 


0.292212075551 


0.023208 


35 


117280195 


0.276622195571 


0.036457 


36 


1142750974 


0.267250650207 


0.019264 


35 


1559953619 


0.285695464272 


0.017605 


34 


3943547917 


0.308882749558 


0.031528 


35 


7346061594 


0.270713039098 


0.023536 


35 


1870287376 


0.278885118559 


0.016743 


33 


8878121745 


0.292567467485 


0.05371 


36 


4559044522 


0.268351164818 


0.016991 


34 


2677976869 


0.290910806654 


0.027865 


35 


0340198507 


0.273008042132 


0.017173 


34 


1860410198 


0.293258906218 


0.029955 


35 


9589514911 


0.273866480307 


0.016559 


34 


3764044862 


0.292368203484 


0.015 


34 


0959952821 


0.145904476604 


0.0544 


36 


9686138698 


0.0771934613908 


0.1561 


39 


294455039 


0.0599663857105 


0.0393 


36 


3055280023 


0.0717497761757 


0.1241 


38 


8100134766 


0.0715971658275 


0.1441 


39 


0528219919 


0.0634435595707 


0.1299 


38 


9956124177 


0.0631822941145 


0.0784 


37 


7810061412 


0.0681153647931 


0.62 


43 


2066240148 


0.374321656282 


0.57 


42 


6299495836 


r\ r\r\ r\ '~7 -i v v r\r\ r\ 

0.329397155099 


0.3 


40 


9068557079 


0.322661461524 


0.38 


41 


981448045 


0.294675002057 


0.43 


42 


2865323661 


0.353944288913 


0.24 


40 


92039446 


0.604691512396 


0.3 


40 


152454556 


0.463180764839 
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0.24 


40.8941528415 


0.281508606621 


0.44 


42.0189541221 


0.290912412354 


0.5 


42.4251183643 


0.285032220696 


0.97 


43.0841865629 


0.742582169563 


0.479 


42.2593258058 


0.305353355969 


0.83 


43.4934896902 


0.372098016098 


0.416 


42.042208973 


0.43961209764 


0.581 


41.9804788609 


0.422581723678 


0.45 


41.7965293235 


0.367900106831 


0.579 


43.1102730617 


0.562403228541 


0.32 


41.1885505108 


0.339690792497 


0.657 


43.0812935967 


0.549444460563 


0.472 


42.0990717933 


0.469994746225 


0.374 


43.1153409984 


0.694417590782 


0.526 


42.013309316 


0.402473462081 


0.763 


44.3657126034 


0.885466302481 


0.58 


42.9398798711 


0.395142768132 


0.43 


41.8899803061 


0.393293335469 


0.45 


42.2825129374 


0.422650328365 


0.828 


44.203018183 


0.611918160502 


0.656 


43.181554143 


0.528708119909 


0.495 


42.1205013485 


0.39827409355 


0.49 


41.8217131258 


0.355251221615 


0.57 


42.742574543 


0.395438890913 


0.388 


42.0969738461 


0.370734294336 


0.45 


42.2128189449 


0.377256221564 


0.48 


42.1266258961 


0.403490417945 


0.615 


42.6108744805 


0.520864721236 


0.4 


42.1077687105 


0.382043588034 


0.655 


42.3305104738 


0.496246176513 


0.498 


43.0296153325 


0.451179870261 


0.465 


41.8375675989 


0.677740680801 


0.453 


42.6855854693 


0.388547543287 


0.425 


41.3010835176 


0.366453610612 


0.278 


40.6729770427 


0.205733686158 


0.477 


42.0071669413 


0.149503540949 


0.95 


43.6972613338 


0.274437861244 


1.057 


44.0787405468 


0.211584422417 


0.816 


43.7951226482 


0.423980625959 


0.455 


42.3522141356 


0.239111591348 
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0.514 


42.4275455763 


0.427580299971 


0.423 


41.6141933692 


0.239548120652 


0.946 


43.1597477534 


0.806741116211 


0.859 


44.022811801 


0.304111603744 


1.031 


42.1761580597 


1.28660173268 


0.936 


44.148887081 


0.518330924179 


0.528 


42.3528772215 


0.240748040733 


0.645 


44.0353749092 


0.964152598465 


0.978 


43.6526492518 


0.312966122482 


0.885 


44.1789952676 


0.466461712419 


0.815 


44.1347089556 


0.919484299622 


0.568 


42.5610152006 


0.293309253351 


0.711 


43.2491247935 


0.248908077637 


0.3396 


41.1022712766 


0.259553392222 


0.3965 


41.4084095732 


0.240552944224 


0.812 


43.8133209255 


0.282936944717 


0.799 


43.1683351205 


0.27444888856 


0.882 


42.9333230617 


0.324847417172 


0.833 


43.407367536 


0.347349605053 


0.874 


43.3865810553 


0.462308366531 


0.772 


43.4354874701 


0.263416996827 


0.543 


42.3919207564 


0.135200439255 


0.75 


43.167266596 


0.135523948437 


0.64 


42.7302552042 


0.225082191179 


0.43 


42.202845059 


0.143459564521 


0.64 


43.1881838551 


0.186543501751 


0.497 


42.4090348089 


0.158488764413 


0.44 


41.9964410659 


0.127761277223 


0.355 


41.516716698 


0.130966220782 


0.78 


43.4907628672 


0.151115530925 


0.54 


42.4339111537 


0.121366799166 


0.86 


43.7451019144 


0.156150591229 


1.02 


44.1993977338 


0.223075604603 


1.14 


44.2839104237 


0.280516457089 


0.854 


43.9153676425 


0.219112489318 


1.37 


44.7910462695 


0.324513218653 


0.975 


44.2284945859 


0.200606912017 


0.97 


44.5842379887 


0.240421078134 


0.74 


43.417393486 


0.196331774539 


1.39 


44.8274132597 


0.591264643225 
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0.46 


42 


1615041182 





190529350474 


1.02 


43 


9810540565 





241216980878 


1.12 


44 


2495126581 





241680223529 


1.23 


44 


852264567 





236250346674 


1.19 


44 


1207308635 





245756009721 


0.839 


43 


5971981016 





211771806558 


1.01 


44 


4466207888 





195202527359 


0.521 


42 


4304300932 





194708817687 


0.475 


42 


1851461163 





216179583718 


0.95 


43 


7698442025 





212959739832 


1.3 


44 


9491044223 





351060245309 


1.305 


44 


4157049301 





292031859582 


0.526 


43 


0211775986 





234516092908 


0.735 


43 


1143549007 





18304597653 


1.14 


44 


0289616697 





770661375912 


1.551 


45 


3317386715 





395907584039 


1.265 


44 


6312640583 





231966862578 


1.34 


44 


953151304 





252440440679 


0.84 


43 


5620294075 





195661579434 


0.468 


42 


4934212383 





156341752732 


0.84 


43 


3513315259 





326964462638 


0.96 


43 


5868914421 





427260522942 


0.8218 


43 


5365504755 





348162814934 


0.93 


44 


6348635492 





689119011672 


0.451 


41 


8469440412 





136544105548 


0.61 


42 


8826976033 





146106297403 


0.83 


44 


6922378489 





457964640542 


0.707 


43 


320290922 





250954794869 


0.415 


41 


7696832434 





12867036751 


0.557 


42 


6475499276 





146252717965 


0.791 


43 


334857994 





19213955201 


0.695 


43 


1559541077 





177154136049 


0.633 


43 


200350196 





204733134997 


0.2486 


40 


6291787524 





144170347426 


0.532 


42 


6510205265 





198544372744 


0.331 


40 


9259072634 





124868930577 


0.346 


41 


2816358992 





133990721793 


0.961 


44 


0860127956 





541232742301 


0.613 


42 


9577096308 





147055302427 


0.3402 


41 


2627973103 





124216129606 


Continued on next page 



88 



Table D.l - continued from previous page 



z 
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0.983 


43.9686072633 


0.873805060546 


0.71 


43.03031003 


0.184976659207 


0.73 


43.3112856025 


0.172061162401 


0.47 


42.2950363729 


0.149449923436 


0.62 


43.0958537208 


0.150521756378 


0.521 


42.2813323021 


0.171703576261 


0.369 


41.5099445593 


0.127128442179 


0.571 


42.4563923477 


0.1847896727 


0.604 


42.4422921595 


0.142594201708 


0.9271 


44.7438233798 


0.616675765778 


0.285 


40.7117195315 


0.128805577715 


0.2912 


40.7936031593 


0.128325808339 


0.548 


42.684557193 


0.163119668307 


0.868 


43.9551353049 


0.697041426681 


0.496 


42.2194365099 


0.160441226732 


0.811 


44.0557629679 


0.356535950064 


0.756 


43.6720421887 


0.2041097647 


0.817 


43.5678420781 


0.305034089394 


0.752 


43.1434994817 


0.179046394457 


0.5516 


42.4072918995 


0.146846411569 


0.3578 


41.3894284101 


0.125749074506 


1.01 


44.6896175735 


0.910412577811 


0.741 


43.4820220235 


0.161255047619 


0.43 


41.8372805299 


0.138494545284 


0.526 


42.6736883867 


0.151797295017 


0.592 


42.524194556 


0.197075800472 


0.905 


43.686095579 


0.424695030791 


0.949 


43.4679668471 


0.467379485606 


0.4607 


41.9741739893 


0.2613454679 


0.3709 


41.5023209499 


0.132033007561 


0.8 


44.4199835555 


0.537991937963 


0.679 


43.4170080051 


0.175601683429 


0.5817 


42.8374976585 


0.184185753799 


0.55 


42.3999711394 


0.158230175273 


0.81 


43.7278491047 


0.491922061811 


0.95 


43.780780395 


0.421021094207 


0.3373 


41.1812895877 


0.124867747916 


0.91 


44.4079361607 


0.394461837584 


0.263 


40.6399477676 


0.123324515918 


0.643 


43.0120301671 


0.161285559543 
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0.691 


43.1407175735 


0.170861512594 


0.357 


41.3588695068 


0.129068537648 


0.721 


43.1285593884 


0.185537803758 


0.581 


42.7918076033 


0.147505707094 


0.6268 


42.7073384153 


0.145364365689 


0.818 


43.5731368588 


0.562299153124 


0.449 


41.9045845939 


0.150334954079 


0.688 


43.0446792667 


0.179014967059 


0.87 


44.0828430243 


0.496674095444 


0.5043 


42.3045856196 


0.141778819311 


0.591 


43.3138714564 


0.267379474385 


0.426 


41.806666569 


0.204450457968 


0.329 


41.3757739861 


0.25327455415 


0.531 


43.0629635087 


0.301071564198 


0.583 


42.5993396847 


0.461245287391 


0.519 


43.0898571036 


0.385243880567 


0.401 


41.6214142647 


0.259723821426 


0.34 


41.1002984666 


0.22330741712 


0.436 


41.8355461879 


0.219396208675 


0.363 


41.3991790454 


0.206870861224 


0.436 


41.9186081822 


0.214365467471 


0.309 


41.2854634774 


0.20604264325 


0.342 


41.3335602868 


0.20174070167 


0.332 


41.2428570494 


0.223960203004 


0.469 


42.0869431125 


0.219485296548 


0.239 


40.4663605358 


0.193168326915 


0.352 


41.4402007529 


0.222632458264 


0.612 


42.2601378946 


0.819156755697 


0.631 


42.5665239435 


0.270098723797 


0.645 


42.8315316986 


0.22803549973 


0.429 


41.8321619796 


0.217936090167 


0.497 


42.0900771628 


0.23056454709 


0.539 


42.3621361424 


0.235072095671 


0.561 


42.5988124145 


0.254763512704 


0.41 


A 'X A c\ r\r\ r\'~y r\ r\r\ 

41.3480927299 


0.217232818458 


0.412 


41.7758022341 


0.2949457048 


0.599 


42.8544760765 


0.290620357209 


0.619 


42.9246968389 


0.246627072238 


0.422 


41.6678602675 


0.238231296443 


0.54 


42.5656285356 


0.254777678043 
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0.401 


42.072813335 


0.312075794951 


0.218 


40.3179435116 


0.198433228858 


0.633 


42.2256486429 


0.296301030677 


0.383 


41.592556266 


0.228094106236 


0.302 


41.5660169963 


0.259288152895 


0.34 


41.1438764437 


0.199964450355 


0.51 


41.9188202391 


0.258309489516 


0.421 


41.965762911 


0.228780937104 


0.399 


41.7087506407 


0.28376407894 


0.493 


42.122916001 


0.229109806841 


0.687 


42.832217021 


0.279368496827 


0.502 


41.8992528194 


0.313490292477 


0.687 


42.8216812976 


0.250279752764 


0.495 


42.2021639825 


0.245006206755 


0.603 


42.6354160597 


0.273892147864 


0.421 


42.0443812883 


0.280749747008 


0.348 


41.5494690313 


0.204875174321 


0.213 


40.3399930846 


0.204646940237 


0.344 


41.0679683623 


0.200210760141 


0.271 


40.5536369049 


0.206465230307 


0.564 


42.3759756991 


0.400947441397 


0.274 


40.7740469723 


0.216212348503 


0.582 


43.2137150031 


0.278992592748 


0.68 


42.8318958632 


0.247902597912 


0.401 


41.9301329836 


0.235368495661 


0.416 


41.7521888571 


0.295167259169 


0.286 


41.387427719 


0.334584228292 


0.314 


41.0209048393 


0.287974294878 


0.581 


43.5994710152 


0.377088903672 


0.463 


41.9546009246 


0.265451577145 


0.341 


41.092049721 


0.209704216095 


0.671 


42.1936025518 


0.331542575911 


0.631 


42.9261055377 


0.232187948068 


0.522 


42.6690298106 


0.304891267733 


0.368 


41.4066563292 


r\ c\r\ A '"T A ^ r\ r\ ^ A r\ r\ 

0.204741881488 


0.309 


40.8507297646 


0.207725830155 


0.528 


42.4787440903 


0.305372522825 


0.268 


40.7695452633 


0.197555926698 


0.695 


42.9207400196 


0.41358331096 


0.284 


40.7849823164 


0.199874464698 
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0.508 
0.781 
0.613 


42.1375516415 
43.3825529133 
42.4302887204 


0.212951542267 
0.341139699918 
0.333098849656 
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Appendix E 
GRB data 



Table E.l: GRB distance measures with la upper and 
lower uncertainties. 



z 


'"p(^) 


a+ 


a 


0.17 


1.0000 






1.036 


0.9416 


0.1688 


0.1710 


1.902 


1.0011 


0.1395 


0.1409 


2.768 


0.9604 


0.1801 


0.1785 


3.634 


1.0598 


0.1907 


0.1882 


4.500 


1.0163 


0.2555 


0.2559 


6.600 


1.0862 


0.3339 


0.3434 
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Appendix F 

Galaxy cluster gas mass fraction data 
set 



Table F.l: Galaly cluster gas mass fraction data set. The 
cluster ID, redshift and the measurement of /gas with 
1(7 error. 



Cluster 


z 




/gas 


± 


(7f 


Abell 1795 





063 





1074 


± 





0075 


Abell 2029 





078 





1117 


± 





0042 


Abell 478 





088 





1211 


± 





0053 


PKS0745-191 





103 





1079 


± 





0124 


Abell 1413 





143 





1082 


± 





0058 


Abell 2204 





152 





1213 


± 





0116 


Abell 383 





188 





0903 


± 





0080 


Abell 963 





206 





1144 


± 





0102 


RXJ0439.0+0521 





208 





0917 


± 





0127 


RXJ1504. 1-0248 





215 





1079 


± 





0111 


Abell 2390 





230 





1257 


± 





0110 


RXJ2129.6+0005 





235 





1299 


± 





0299 


Abell 1835 





252 





1197 


± 





0082 


Abell 611 





288 





1020 


± 





0133 


Zwicky 3146 





291 





0943 


± 





0163 


Abell 2537 





295 





0949 


± 





0147 


MS2137.3-2353 





313 





1106 


± 





0061 


MACSJ0242.6-2132 





314 





1268 


± 





0131 


MACSJ1427.6-2521 





318 





1052 


± 





0220 


MACSJ2229.8-2756 





324 





1452 


± 





0265 


MACSJ0947.2+7623 





345 





1048 


± 





0196 


MACSJ1931.8-2635 





352 





1193 


± 





0266 
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Cluster 


z 


/gas ± lO"! 


MACSJ1115.8+0129 





355 





0925 


± 





0283 


MACSJ1532.9+3021 





363 





1280 


± 





0162 


MACSJOOl 1.7-1523 





378 





1067 


± 





0125 


MACSJ1720.3+3536 





391 





1153 


± 





0151 


MACSJ0429.6-0253 





399 





1375 


± 





0154 


MACSJ0159.8-0849 





404 





1097 


± 





0160 


MACS J2046. 0-3430 





423 





1253 


± 





0398 


MACSJ1359.2-1929 





447 





0845 


± 





0290 


MACSJ0329.7-0212 





450 





1262 


± 





0129 


RXJ1347.5-1144 





451 





0923 


± 





0078 


3C295 





461 





1067 


± 





0096 


MACSJ1621.6+3810 





461 
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